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The  Reattachment  of  Laminar  Cavity  Flow  with  Heat  Transfer 


At  Hypersonic  Speed 
Summary 

A  study  of  cavity  flow  is  made  based  upon  the  analysis  of  a  physical  model 
and  by  referring  to  the  experimental  data  of  a  particular  test  model  with  the  reat¬ 
tachment  surface  perpendicular  to  free  stream  direction  at  free  stream  Mach 
number  of  eleven. 

The  pressure  and  heat  transfer  in  the  reattachment  zone  may  be  predicted 
approximately  by  the  flow  quantities  along  the  dividing  stream  line  employing  the 
mixing  theory  and  stagnation  heat  transfer  equation  for  a  blunt  body.  Good  agree¬ 
ment  with  experimental  data  is  obtained. 

It  is  found  that  the  value  of  the  average  mixing  rate  correlation  function  for 
the  cavity  flow  is  15  and  is  equal  to  that  value  determined  by  Glick  for  separated 
flow  caused  by  shock  impingement.  Semi-empirically,  efficiency  factor  of  com¬ 
pression  equal  to  0.  67  at  the  reattachment  is  evaluated  and  the  reattachment 
velocity  gradient  parameter  for  the  particular  test  model  is  similar  to  this  case  of 
stagnation  flow  over  blunt  nose.  It  appears  that  the  heat  transfer  at  the  reattach¬ 
ment  surface  perpendicular  to  the  free  stream  direction  may  be  predicted  approx¬ 
imately  by  the  stagnation  point  heat  transfer  equation  for  a  two-dimensional  blunt 


surface. 
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Introduction 

The  problem  of  cavity  flow  which  involve  s  separation  as  well  as  reattac' 
ment  is  complicated  and,  as  yet,  not  completely  understood.  The  solution  of 
this  problem  is  not  only  for  the  interest  of  basic  science  but  also  important  fo 
many  practical  applications.  For  example,  for  a  cut-out  in  a  structure  surfa> 
such  as  a  bomb  bay,  open  cock  pit,  escape  hatch,  etc.  ;  the  heat  transfer  at  th 
reattachment  zone  may  become  very  high  at  the  hypersonic  flight  speeds. 

The  pronounced  features  of  cavity  flow  are  viscous  momentum  transfer 
along  the  boundary  stream  line  of  the  cavity  opening  and  the  circulating  flow 
within  the  cavity.  The  momentum  is  transferred  from  the  external  flow  to  the 
internal  flow  due  to  viscous  mixing  and  the  flow  within  the  cavity  circulating  t( 
balance  the  mass  rate  flow. 

For  subsonic  flow  the  concept  of  mass  balance  between  the  outgoing  and 
incoming  fluid  from  and  into  the  idealized  wake  model  has  been  used  by  Tanne 
to  predict  with  a  simple  analysis  the  drag  and  base  pressure  of  cones  as  well  ; 
the  area  of  outgoing  flow  from  the  wake.  However,  no  consideration  is  given  i 
the  viscous  momentum  transfer. 

For  supersonic  laminar  separated  and  reattached  flow  Chapman  [2,  3]  « 
Crocco-Lees  [4]  developed  the  detailed  mixing  theory  for  the  interaction  betw 
viscous  internal  flow  and  the  nearly  isentropic  external  stream.  Chapman  [2 
troduced  the  idea  of  the  so-called  "dividing  stream  line"  which  divides  the  flov 
field  within  the  very  deep  cavity  into  two  zones,  i.  e.  ,  outside  by  pass  flow  zon 
of  this  stream  line  and  inside  zone  of  circulating  flow,  and  obtained  an  exact 
similarity  solution  for  the  velocity  profile  <  f  the  mixing  layer.  A  simplified 


calculation  showed  that  the  reattachment  process  is  essentially  isentropic. 

The  Crocco-Lees'  theory  is  general  and  applicable  qualitatively  for  a  wide 

range  of  separated  flow  problems.  The  separated  and  reattached  flows  are 

characterized  by  two  governing  parameters  of  <  and  f  defined  by  boundary 

layer  thickness.  Much  of  the  complex  phenomena  of  the  separated  and  reattached 

flow  are  clarified  by  this  theory.  Glick  [5]  used  Chapman's  flow  model  which 

incorporates  the  concept  of  the  "dividing  stream  line"  to  apply  the  Crocco-Lees' 

theory.  This  model  has  been  translated  into  Crocco-Lees  language  by  the  semi- 

empirical  approach.  Glick  found  that  the  cause  of  poor  quantitative  prediction  of 

Crocco-Lees'  theory  was  an  incorrect  mixing  rate  correlation  function  C  (<), 

thus  by  determining  the  proper  approximate  C  (<)  -  F  (  K)  relation  where  F  =  ~Z  ~ 

< 

for  separated  and  reattached  flow  a  good  agreement  with  experimental  data  was 
obtained  for  the  pressure  distribution  in  the  region  of  separated  flow  caused  by 
shock  interaction. 

Restricting  to  a  moderately  compressible  separated  flow  over  the  concave 
surface.  Bloom  [6]  obtained  solutions  for  velocity  gradient,  pressure,  boundary 
layer  thickness  employing  an  integral  method  and  emphasizing  the  reverse  flow 
profiles. 

For  the  supersonic  turbulent  reattachment  flow  behind  a  step,  Chapman- 
Korst  [7,  8]  established  a  recompression  criterion  of  reattachment  and 
McDonald  [9]  improved  this  criterion.  McDonald  formulated  a  transformed  shape 
parameter  for  the  boundary  layer  by  using  Stewartson-Illingworth  transformation 
to  specify  the  reattachment. 

For  the  solution  of  slow  circulating  flow  within  the  cavity,  Squire  [10  ] 
proposed  an  analytical  approach  to  divide  the  cavity  zone  into  two  parts;  the 
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inside  part  of  "core"  and  the  outside  part  of  the  boundary  layer  which  surrounds 
the  core.  The  motion  in  the  cavity  is  maintained  by  the  shear  stress  of  the 
outer  flow  acting  on  the  stream  line  boundary  of  the  cavity.  Mills  [  11  ]  obtainec 
analytically  the  velocity  profiles  along  the  square  cavity  wall  based  upon  Squire 
idea  and  his  results  are  in  a  good  agreement  with  the  experimental  data. 

Batchelor  [12]  proposed  a  model  of  a  finite  laminar  separated  flow  regioj 
for  a  two-dimensional  as  well  as  axial  symmetric  cavity  with  limiting  steady 
flow  at  high  Reynolds  number.  In  this  model  an  inviscid  rotational  core  is  sep¬ 
arated  from  the  external  flow  by  a  recirculating  viscous  shear  layer.  Since  thi 
cavity  model  does  not  provide  the  unique  flow  field,  the  recirculating  flow  is 
divided  further  into  a  cellular  substructure  of  eddies  each  satisfying  the 
Batchelor's  conditions. 

Burggraf  [13]  applied  this  Batchelor's  model  to  analyze  the  flow  and  heal 
transfer  in  rectangular  cavities  of  various  depths  by  employing  an  approximate 
linear  theory  in  order  to  account  for  the  effect  of  cavity  depth  by  a  correlation 
factor.  For  convenience,  the  problem  is  broken  up  into  several  parts  and  eacl 
is  analyzed  separately,  then  by  matching  these  partial  solutions  the  complete 
solution  is  synthesized.  Burggraf's  [13  ]  incompressible  flow  analysis  can  be 
extended  to  compressible  flow,  but  his  analysis  predicts  a  pressure  lower  thar 
that  obtained  by  experiment  at  the  reattachment. 

Most  recently  Burggraf  [14  ]  car  ried  out  an  extensive  analysis  of  incom¬ 
pressible  flow  for  the  viscous  structure  of  a  separated  eddy  within  a  fixed  finit 
cavity  as  a  function  of  the  Reynolds  number.  In  particular  he  considered  the 
square  cavity  and  determined  the  velocity  profile,  distributions  of  skin  frictioi 
and  total  power.  The  thermal  condition  in  a  recirculating  eddy  is  considered, 


but  no  specific  information  is  available  for  the  pressure  and  heat  transfer  at  the 
reattachment. 

On  the  heat  transfer  of  the  laminar  cavity  flow  Chapman  [3  ]  succeeded 
only  to  compute  the  average  heat  transfer  rate  amounting  to  0.  56  of  that  of  the 
attached  flow.  For  the  determination  of  the  local  heat  transfer  rate  the  flow 
solution  which  considers  the  effects  of  the  circulating  flow  and  the  relation  of 
flow  direction  with  respect  to  surface  configuration  should  be  known.  Chung  and 
Viegas  [15  ]  determined  analytically  the  velocity  distribution  along  the  reattach¬ 
ment  wall  of  the  cavity,  assuming  that  the  reattachment  flow  is  laminar,  two- 
dimensional  as  well  as  incompressible  but  rotational.  They  further  assumed 
that  the  reattachment  of  the  flow  along  the  dividing  stream  line  takes  place  normal 
to  the  reattachment  surface.  By  using  their  analytical  results  and  a  semi-empirical 
formula  for  heat  transfer,  an  average  value  of  the  heat  transfer  in  the  reattachment 
zone  was  obtained.  Burggraf’s  [13]  analytical  solution  at  the  reattachment  results 
in  the  square  root  singularity. 

Carlson  [  16  jdetermined  the  local  heat  transfer  for  laminar  separated 
flow  upstream  of  reattachment  behind  a  step  of  a  body  of  revolution  at  hypersonic 
speed  taking  into  account  a  velocity  profile  of  reverse  flow.  It  has  been  assumed 
that  no  chemical  reaction  and  no  mass  addition  of  fluid  takes  place,  the  cavity 
pressure  is  constant  throughout  the  region,  and  the  dividing  stream  line  is  straight. 
Although  it  was  shown  that  the  rate  of  heat  transfer  increases  approaching  reattach¬ 
ment,  no  detailed  information  on  the  reattachment  heat  transfer  is  available. 
Recently  Korst  [  1  7 ]  developed  the  generalized  analysis  of  dynamics  and  thermo¬ 
dynamics  of  separated  flow. 
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This  review  of  references  shows  that  up  to  date  no  reliable  analysis  for 
the  laminar  local  properties  of  reattachment  flow  and  the  heat  transfer  for  the 
cavity  has  been  developed.  Therefore,  for  the  first  step  of  the  investigation 
based  on  two  flow  models  of  Chapman  and  Squire  an  attempt  is  made  to  find  out 
whether  the  existing  analysis  can  be  applied  for  reattachment  of  the  cavity  flow. 

Various  experiments  on  the  flow  and  heat  transfer  of  cavity  have  been 
carried  out  in  the  past.  In  the  region  of  subsonic  speeds,  Roshko  C  1 8 3  measure 
pressure  and  velocity  distributions  on  the  wall  of  rectangular  cavity  by  varying 
the  ratios  of  depth  to  length  of  cavity.  Tani  [  1 9  ]measured  longitudinal  compon 
of  mean  velocity,  turbulent  intensity  and  turbulent  shear  stress  across  several 
traverse  sections  of  rectangular  cavity.  Maull  and  East  [20]  discovered  exper: 
mentally  over  the  bottom  of  rectangular  cavities,  the  three-dimensional  cellula: 
flow  pattern  between  certain  range  of  ratio  of  cavity  length  to  width,  but  Fox  [2 
did  not  notice  it.  However,  Fox  [21]  indicates  that  there  exists  the  distinct  flo 
regime  in  the  pressure  coefficient  and  velocity  profile  and  used  them  to  identify 
the  range  of  ratio  of  length  to  width  of  cavity  with  different  flow  regime. 

In  the  range  of  supersonic  speeds,  by  using  rectangular  cavities,  Charw 
[22]  evaluated  the  pressure  distribution  on  the  floor  and  recompression  surface 
as  well  as  velocity  distribution  through  shear  layer.  Furthermore,  the  variatic 
of  heat  transfer  distribution  on  the  floor  was  measured.  Larson  [23]  determin< 
for  laminar  and  turbulent  flow  the  average  heat  transfer  for  the  cavity  formed 
on  an  axially  symmetric  body  at  free  stream  Mach  numbers  up  to  4  and  also  the 
streamwise  distribution  of  heat  transfer  including  the  reattachment  zone. 
Bogdonoff  and  Vas  [24]  measured  distributions  of  pressure  and  heat  transfer 
over  the  cone-cavity  model  at  M^  ^  11.  7.  More  recently,  Nicoll  [  25]  carried 
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out  extensive  measurements  on  streamwise  pressure  and  heat  transfer  distri¬ 
butions  for  cut-outs  on  cones  using  helium  as  a  fluid  medium  at  free  Mach 
numbers  of  11  and  20.  These  results  show  that  Chung-Viega°'  results  of 
average  reattachment  heat  transfer  are  about  twice  those  obtained  in  these  ex¬ 
periments.  The  existing  experimental  data  are  useful  to  evaluate  semi- empirically 
reattachment  flow  parameters  such  as  reattachment  velocity  gradient  and  efficiency 
factor  of  compression  and  confirm  the  applicability  of  the  existing  analysis  to 
cavity  flow. 
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1.  Physical  Model  of  Cavity  Flow 

For  the  investigation  of  flow  in  the  region  of  the  cavity,  a  physical  model 
as  shown  in  Figure  1  is  proposed.  This  model  divides  the  whole  flow  field  into 
two  regions  by  the  dividing  stream  line,  i.  e.  ,  the  external  nearly  isentropic 
flow  region  and  the  internal  viscous  flow  region.  This  represents  basically  th< 
Chapman’s  model  [2  ]  as  sketched  in  Figure  1  except  for  the  shaded  area  withii 
the  cavity.  This  model  emphasizes  the  viscous  mixing  process  of  the  external 
flow  with  the  internal  flow. 

The  velocity  profile  at  the  separation  point  S  is  similar  to  that  at  the  rcat- 
tachment  point  R  where  shear  stress  is  zero.  The  fluid  particle  which  is  adja 
cent  to  the  wall  at  the  separation  must  be  adjacent  to  the  wall  at  the  reattachm 
and  the  external  flow  may  be  considered  as  an  essentially  by-pass  flow.  Due  t 
the  viscosity  of  the  fluid,  mixing  of  external  fluid  with  internal  fluid  takes  plac 
in  the  mixing  layer  which  extends  from  the  separation  point  to  the  beginning  of 
reattachment  zone.  But  with  the  presence  of  heat  transfer,  at  supersonic  spe< 
the  process  of  mass  exchange  may  be  considered  as  quasi-steady,  and  the  cor 
responding  values  at  the  proper  flow  properties  are  taken  for  the  analysis,  si 
the  dividing  stream  line  may  fluctuate  and  periodic  mass  transfer  can  take  pla 
as  evidenced  by  Charwat  [22]  for  the  rectangular  cavity.  The  thickness  of  the 
thin  mixing  layer  grows  parabolically  with  downstream  distance  from  the  orig 
of  mixing  but  the  rate  of  growth  is  about  three  times  larger  than  that  of  the 
attached  laminar  layer. 

The  internal  flow  field  is  a  circulating  flow  region  and  this  flow  field  may 
be  divided  again  into  two  regions,  i.  e.  ,  the  inner  "core"  and  the  outer  thin 
boundary  layer  belt  surrounding  the  core  as  proposed  by  Squire  [10]  for  subs< 
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speed.  Figure  2  illustrates  the  recirculating  cavity  flow  while  the  shaded  area 
of  Figure  1  indicates  the  core.  Within  the  core  the  mass  entrained  from  the  so- 
called  dead  air  region  where  the  velocity  is  small  but  not  necessarily  zero,  (i.  e.  , 
the  region  below  the  zero  velocity  line)  and  the  mass  reversed  back  into  the  dead 
air  region  to  balance  each  other  and  are  responsible  for  the  maintenance  of  the 
circulatory  flow  within  the  cavity.  The  back  flow  is  caused  by  the  pressure  rise 
in  the  reattachment  zone  [Figure  3]  and  the  mass  balance  causes  the  pressure 
of  cavity  to  differ  from  that  upstream  of  separation. 

Recent  a.ialytical  investigation  of  Denison  and  Baum  [26]  by  solving  numer¬ 
ically  the  governing  differential  equations  on  a  computer  showed  that  sufficiently 

if 

downstream  of  the  separated  flow  region  u  =  ud'ue*  defined  as  the  ratio  of 
laminar  flow  velocity  along  the  dividing  stream  line  and  the  potential  stream  line 
outside  the  cavity,  is  u  =  0.  587  regardless  of  whether  it  is  axi-symmetric  or 
two-dimensional  flow.  This  result  is  obtained  by  considering  an  initial  boundary 
layer  thickness  at  the  separation.  On  the  other  hand  for  two-dimensional  laminar 
flow.  Chapman  [2,7]  obtained  analytically  the  same  value  of  u*  under  the  assump¬ 
tion  of  zero  boundary  layer  thickness  at  the  separation.  Hence,  it  may  be  con¬ 
cluded  that  if  the  reattachment  distance  is  small  compared  to  the  whole  cavity 

length,  the  reattachment  process  can  be  investigated  separately  apart  from  the 

i 

upstream  conditions  where^the  momentum  transfer  up  to  the  initial  pressure  rise 
for  reattachment  is  known. 

In  the  reattachment  zone  the  pressure  rises  as  shown  in  Figure  3  because 
ot  the  increased  momentum  in  the  shear  layer  due  to  viscous  mixing,  and  the 
deceleration  of  the  flow.  Since  at  the  reattachment  the  velocity  is  zero,  a 
particle  moving  along  the  stream  line  within  the  mixing  layer  must  be  able  to 
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overcome  the  pressure  rise  and  therefore  in  order  to  pass  through  the  down¬ 
stream  direction,  its  total  pressure  at  the  reattachment  must  be  larger  than 
p  the  static  pressure  downstream  of  the  reattachment.  Thus  although  partic! 
(a)  (see  Figure  3)  may  pass  through  reattachment,  particle  (b)  may  not,  becauj 
its  velocity  is  lower  than  that  corresponding  to  the  pressure  at  the  reattachmer 
equal  to  p  .  Therefore,  particle  (b)  is  reversed  back.  Chapman  et  al  [7  ] 
observed  that  the  reattachment  pressure  rise  is  independent  of  Reynolds  num¬ 
ber.  Furthermore,  as  indicated  in  Figure  4,  Harper  [27]  found  that  if  the  in¬ 
compressible  flow  stagnates  on  a  flat  surface  then  the  flow  region  in  the  neigh¬ 
borhood  of  the  stagnation  point  is  inviscid,  although  the  stagnation  point  itself 
lies  in  the  viscous  layer.  Hence,  the  reattachment  process  approaching  the 
reattachment  point  may  be  considered  as  isentropic  but  in  the  very  close  regioi 
near  the  reattachment  as  well  as  at  the  reattachment  point  itself  viscosity  effec 
should  be  included.  Besides  the  effect  of  compressibility,  geometrical  relatioi 
for  the  reattaching  stream  line  to  the  reattaching  wall,  etc.  are  also  to  be  con¬ 
sidered  to  evaluate  the  flow  properties  at  the  reattachment  correctly. 

Similar  to  the  stagnating  flow  process,  the  flow  velocity  reduces  to  zero  al 
the  reattachment  point  but  its  velocity  gradient  is  finite.  The  flow  velocity  aloi 
the  reattaching  surface  may  be  computed  based  upon  the  cavity  flow  model  in 
the  core  and  surrounding  boundary  layer.  Since  heat  transfer  in  the  stagnation 
region,  can  be  evaluated  by  knowing  the  velocity  gradient  at  the  stagnation  poin 
and  the  velocity  distribution  around  it,  the  heat  transfer,  in  the  reattachment 
zone  may  also  be  solved  based  upon  the  flow  model  proposed  here. 
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2.  Pressure  Rise  in  the  Reattachment  Zone 

The  pressure  rise  in  the  reattachment  zone  of  the  cavity  is  determined 
analytically  and  its  numerical  values  are  evaluated  semi-empirically. 

The  reattachment  pressure  rise  is  due  to  momentum  transfer  upstream 
of  reattachment  and  deceleration  in  the  reattachment  zone.  The  mass  transfer 
is  computed  by  the  method  of  Glick  [5]  and  the  deceleration  of  flow  by  Harper's 
[27]  formula. 

The  coordinate  system  is  indicated  in  Figure  5.  The  origin  of  the  coordi¬ 
nate  is  located  at  the  separation  point  with  the  x-axis  along  the  dividing  stream 
line  and  the  y  axis  perpendicular  to  x-axis.  If  the  coordinate  along  the  zero 
velocity  line  [Figure  1]  is  designated  by  x,  which  corresponds  to  imaginary 
wall-coordinate  then  the  direction  of  x  is  approximately  equal  to  that  of  x  down¬ 
stream  of  the  separation  point  in  the  mixing  zone.  The  initial  point  of  the  pres¬ 
sure  rise  for  the  reattachment  is  denoted  by  subscript  o. 

For  the  analysis  the  following  assumptions  are  made: 

1.  Cavity  is  open,  i.e.,  the  dividing  streamline  bridges  the  cavity  from 
separation  to  reattachment  and  does  not  attach  on  the  cavity  floor.  L/D 
is  sufficiently  small  so  that  the  streamwise  pressure  gradient  from 
separation  to  the  initial  point  of  the  pressure  rise  is  zero.  The  location 
of  the  initial  point  for  reattachment  pressure  rise  is  considered  known. 

2.  Pressure  gradient  traverse  to  the  stream  direction  is  zero. 

3.  Flow  is  quasi-steady  and  gas  is  thermally  as  well  as  calorically  perfect. 
The  property  values  of  flow  are  interpreted  properly  as  such. 

4.  The  gradient  of  viscious  or  Reynolds  stress  in  the  flow  direction  are 
negligible  compared  to  the  static  pressure  gradient  in  the  flow  direction. 
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5.  Theory  of  Crocco-Lees  for  two-dimensional  flow  is  applicable  for 
axially  symmetric  flow. 

6.  Heat  transfer  effect  is  negligible  for  the  computation  of  pressure  rise 

7.  The  pressure  rise  prior  to  the  initial  point  of  reattachment  is  negligit 

For  a  streamline  within  the  mixing  layer  with  angle,  _  @  t  between 

Ox 

streamline  and  the  outer  boundary  of  the  dissipative  flow,  at  y  =rf ,  [Figure  ' 
the  governing  equations  are; 

dm  r  d/  .  a  .  ,  , 

dx  =  e  ue  L  cfc  0  ]  =  k  .  e  ue  (j) 

where  k  =  d  &  -  6  is  mixing  coefficient  and  0  is  streamline  direction  anj 

dx 

relative  to  the  wall  at  y  =  J 

This  equation  expresses  the  rate  at  which  mass  is  transferred  from  exte 
flow  to  internal  flow  due  to  mixing  and 


where 


ue 


dm  ,  dp 

d  x  )  ’  d  dx  Xw  (2) 

--  momentum  flux  in  x-direction 
--  mass  flux  in  x-direction 


o 

is  the  momentum  equation  in  x-direction. 

These  definitions  are  referred  to  the  limit  o  and  <f  ,  i.e.,  the  partially 
truncated  thickness  of  the  viscous  layer  measured  from  the  axis  to  the  dividir 
steamline  of  the  outer  edge  of  the  viscous  layer  as  shown  in  Sketch  A 

By  r  eferring  to  the  total  thickness  of  the  viscous  layer  cTt,  as  shown  in 
Sketch  B,  it  may  be  postulated  that  the  mass  flow  rates  for  both  (profiles  as 
shown  in  Sketch  A  and  Sketch  B)  remain  the  same,  although  J <  c/^. 

Consider  two  different  portions  of  the  flow  path.  The  first  consists  of  a 
constant  mass  rate  of  circulatory  flow  along  b  path  below  the  zero  velocity  lir 
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(Sketch  C),  say  mc  and  another  consist  of  the  mass  rate  of  the  by-pass  flow 
over  the  zero  velocity  line,  say  mp  which  corresponds  to  that  indicated  by 
Sketch  B . 

Now  consider  the  divided  flow  fields  above  and  below  the  dividing  streamline. 
Then  due  to  the  viscour  mixing  effect  between  the  dividing  stream  line  and  zero 
velocity  line  and  the  breathing  of  cavity,  an  additional  mass  Am  is  entrained  and 
this  am  flows  out  from  the  cavity,  along  path  a  (Sketch  C)  in  order  to  maintain 
mc  inside  the  cavity.  Hence,  am  is  ejected  from  the  cavity  into  the  viscous  mix¬ 
ing  layer  as  observed  by  Charwat  et  al  [22]  for  the  supersonic  ta .ity  flow.  Tani 

[2  8]  observed  also  the  mass  ejection  from  the  separated  flow  region  behind  a  step 
across  the  dividing  sti  earn  line  at  subsonic  free  stream  velocity.  Hence  at  the 

same  distance  from  the  top  of  profile  of  Sketch  A  and  B  in  negative  y  direction  it 
maybe  assumed  that  P.u  >  p.'u'andthe  mass  flow  rate  over  the  dividing  stream 

r* 

line  in  the  region  of  shortened  viscous  layer  j  f*  u  dy  as  shown  in  Sketch  A  is 

r"t  o 

postulated  to  be  equal  to  I  ^'u'dy  which  is  based  on  a  hypothetical  profile 

o 

with  no  entrainment  from  the  inner  cavity  where  is  its  actual  thickness 

of  the  viscous  layer. 

The  definition  of  momentum  flux  in  x-dirertion  I  is  made  by  an  integration 
from  dividing  stream  line  to  </  ,  similar  to  the  mass  flux  in  x-direction  m  . 

For  the  external  isentropic  flow,  the  Bernoulli's  equation  is, 


where 


I  !e 

p  dx 

Y  e  -  f 


1  dwg 

Te  dx 


(3) 


y  ~  cp 

Wt  =  ue  /  at 


1  -  ( 

/  c  v 


Y  -  1 


we 


} 


2)1  /  Y 


and 


=  /7  X’  T, 


stagnation  speed  of  sound. 


In  order  to  compute  the  momentum  transfer  in  a  region  between  separat: 


and  the  beginning  of  reattachment  and  pressure  rise  in  the  reattachment  zone 
these  equations  are  cast  into  the  language  of  Crocco-Lees  [4]  using  Glick's  [ 


improved  method. 


A  basic  non-dimensional  shape  parameter  is  defined  by 

I  momentum  flux 

K  *  IE  *  mass  flux  times  local  external  velocity 

m  Ue 


or  <  = 


actual  momentum  flux 

momentum  flux  of  mass  flux  moving  at  u  =  ue 


where  subscript  1  refers  to  average  value  of  viscous  region. 

This  parameter  <  has  the  ability  to  correlate  the  velocity  profile  in  the  pi 
sure  plateau  region  of  the  separated  flow  which  is  different  from  the  known  Bis 


flow. 


The  two  key  correlations  of  F  (  <)  and  C  (<  )  for  the  solution  of  the  separai 


and  reattached  flow  are  formulated  by  <  as  follows: 


F  (  <  )  =  *  -  1 

X  & 


c  (  <  )  = 


k  m 


where  f  =  1 1  + 

Tt 

\  /  /  0 „  uot 


a  2 
1C  w 


(dm/dx)  ’  (jeue) 
)i  e/  m 


T,  §  ,  n  are  temperature,  density  and  dynamic  viscosity  respectively  and 
the  subscript  t  refers  to  total  temperature. 

Then  the  continuity,  momentum  and  Bernoulli's  equations  are  reduced  to 


P 

(  m  k  we)  =  wE 
d  we 


_  P  we 


where 


m  =  m  •  at 
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Since  upstream  of  reattachment,  ^EL  and  cf  are  negligible,  momentum 
equation  (5)  reduces  to 


(  m  k  we)  - 


1  -  <  dx 


J_  dm 

m  dx 


By  integrating  from  separation  point  to  the  initial  point  of  reattachment  zone. 


<  =  1  -  (  1 
o 


K  \  mB 


where  the  subscript  s  and  o  refer  to  separation  and  initial  point  of  reattach¬ 
ment  respectively.  From  this  equation  it  is  seen  that  the  levels  of  momentum 
and  energy  in  the  viscous  region  are  raised  by  mass  entrainment,  because  for 

mo  >  ms »  <  o  >  1 . 

Since 


dm  _  p  ^  _  C  ( K  ) 

% - » 

and 

JB_  =  P  = 

'  Tj  _1_  P  Y 

Tt  y  w  e 

it  follows 


IL  u 
f.  "eat 


At  u 

P  Y  -=r  .1 

Te  »t 


_  e  e  t 


dm  _  C(<  ) 

dx  m“ 


? e  ue  Me  at 


The  increase  of  mass  flow  due  to  entrainment  can  be  evaluated  by  integrat¬ 
ing  this  equation,  if  the  relation  of  C  (<  )  with  respect  to  x  is  known.  Although 
it  is  expected  that  C  (<)  starts  at  zero  and  rises  near  separation,  for  simplicity 
an  average  constant  value  of  C  is  used  and  its  numerical  value  is  determined 
properly  from  experimental  data. 


Then  by  integrating  equation  (8), 


m_  -  m . 

o  Ss 


2  C  fe  ue  Pe  at  *o 


2 


or 


3 

™s 


7  = 


2  R, 


-  2  2 
>xo  C  at  Me 
~L 
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(9) 


8 


Where  xQ  is  streamwi9e  distance  from  separation  to  the  beginning  of  reattach¬ 
ment  and  Rexo  =  f  eUeX° 


Taking 


m 


s  * 


Re 


Mt 


1  -  < 


and  Re6**  =  \J  A  \J  Re} 


1  - 


x 

1  -  * 


(Ref.  30) 
(Ref.  5) 


/TT 


2(1  -<»)*  c  .  X£> 
A  ^ 


The  value  of  A  may  be  taken  A  =  0.  44  and  l  is  the  distance  from  the  sta 
nation  point  to  the  separation  point. 

Since  in  the  reattachment  zone  the  streamwise  pressure  rise  may  be  com¬ 
puted  by  flow  quantities  along  the  dividing  stream  line  [7],  and  mixing  rate  and 
viscous  term  are  negligible,  the  momentum  equation  (5)  is  reduced  to 


t  -  -  6  g 

where  the  subscript  d  refers  to  dividing  stream  line 
Since  m  =  p  6/  , 


(10) 


it  follows 


*1  K  +  K 


d  w. 


MV. 


W  d 


d£ 

P 


but 


thus 


fi  - 


V  Wd 


(  F  +  Td 


_?  ) 

•t 


dw, 


K  —  - 


*  - ^  +  K  [  F 

W.  L 


1  ] 


dw , 


MV, 


Because 


dwd 


(  Tt  )  £i*d 


Td 


Md 


finally 


<1  <  =  < 


dMc 

m7 


(11) 


■*» 


T 
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In  order  to  integrate  this  equation,  a  relation  of  F  with  respect  to  *  must 
be  known. 

In  the  reattachment  zone,  a  linear  relation  of  F  with  respect  to  <  may 
be  assumed  [5], 
thus 

F  ■  or  k  +  g  (12) 

where  a  and  g  are  constants  depending  on  the  values  of  F0  and  *  Q  . 
Hence  equation  (11)  becomes 


_ A<  ,  _  =  d  Md 

Integrating  from  the  beginning  point  of  reattachment  to  downstream  distance. 

r  <  _  -  I 

Md0  =  Md  I  _JL  (  «  +3  /  a  *  j  J 

L  ‘  <  o  +  J  /  a  J  (13) 

Putting 

0  =  F  /  K 

rj  =  Fq/F 

«)  =  «£,/  K 
it  follows ; 

=  Q(rj-l)  /  (  a>  —  1 ) 

g  =  F  (u>-r|)/(uj-l) 

and  equation  (13)  is  rewritten  by 

-  (*  -  1) 

F  (r,  -  a)  (14) 

«do  ■  Md  (tj) 

The  value  of  FQ  is  taken  as  F0  =  Fs  =  2.  85  [5]  . 

Assuming  the  isentropic  compression,  the  pressure  rise  in  the  reattachment 


zone  is  given  by 


cav 


V  -  1 

1  +  -y—  Md  (  =  , 


2  (D 


-2  (a;  -  1) 
F(tj-  u,) 


-1 


(15) 


1  +  V  m  l-  Md 

In  order  to  evaluate  the  local  pressure  rise,  a  trajectory  of  F  -  <  such 
as  shown  in  Figure  7  is  drawn.  For  completion  of  this  trajectory,  since  the 
value  of  F  and  <  at  separation  and  reattachment  points  are  known,  only  those 
values  of  F  and  <  at  the  beginning  of  reattachment  should  be^omj^ijted.  Usini 
the  completed  F  -  K  trajectory  local  values  of  ^  are 

evaluated  at  the  proper  locations  in  the  reattachment  zone. 

The  variation  of  Md  in  the  reattachment  zone  is  determined  by  using  Harp 
[2  7]  analytical  solutions  of  velocity  distribution  approaching  the  stagnation  poi 
The  velocity  of  flow  approaching  the  stagnation  point  depends  on  vortices  |"~'0 

r 1  of  inside  and  outside  the  stagnating  stream  lines  and  angle  0  as  indica 
in  equation  (16).  Harper's  analysis  is  applicable  for  incompressible  two-dim 
sional  flow  but  it  is  assumed  that  his  solution  is  also  applicable  for  axi-symmi 

reattaching  flow. 


Harper's  solutions  are: 


for 

6  <  T 

*  > 

ud  (x)  =  2  ( 

1  o  +il) 

(L  -  x)  tan  0  +  0  (L  -  x) 

TT 

0  =  j 

ud  (x)  =  -  (  / 

-  ,  r 

o  +  /  1) 

"  ]2  L  '  a]  ’  +  ‘j  <16> 

+  0  (L  -  x) 

and 

3  >  j 

ud  (x)  =  H 

0 

C  (L  -  x) 

n  ”  1 

^  +  j  <ro+ri)  (L  -  x)  tan  0 

where  C  and  a] 

are  constants. 

+  0  (L  -  x) 

It  is  noticed  that  for  0  -  II  ,  the  inviscid  solution  of  stagnating  flow  has 
a  logarithmic  singularity. 


Ma  is  evaluated  by  Busemann's  integral  [30]  for  a  perfect  gas,  taking  the 
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cavity  gas  temperature  equal  to  recovery  temperature. 


Since 


Ad  At 

17  Te 


Y  -  1 
2 


=  1  + 


-)2  { 
ue  1  ; 


and  Td 

t; 

it  follows 


V  -  1 


1  +Y_ 


Md2 


[e  (  ud/ue)< 


«,•] 


By  denoting 


ud  _ 


2  *2 


u  /  ( 


1  *  Y-4-i  Me2  (1  -  u*2)  j 


Then  the  pressure  at  the  reattachment  point  is 


■  f  1  ■  *  v 


Mo  ) 


Y 

r- 1 


1  +  (1  -  u*2)  Y 


Y 

—  1 


J  is  an  empirical  efficiency  factor  of  compression  relative  to  that  of  isentropic 

1 

process  and  J  <  1.  If  ^  =  1,  then  compression  takes  place  isentrop- 


ically. 


Numerical  Evaluation  of  Pressure  Rise  in  the  Reattachment 


Zone 


The  value  of  ,  is  evaluated  based  on  the  experimental  data  of  a  parti 

ular  test  model  of  Nicoll  [25]  as  shown  in  Figure  8.  With  this  model  the  ini 

point  of  rise  of  pressure  and  heat  transfer  is  located  approx imately  at  the  s< 

distance  of  1  -  =0.2.  By  trial  and  error,  the  value  of  C  is  determinet 

Jb 

which  yields  the  pressure  rise  in  the  reattachment  to  be  in  good  agreement  \ 

experimental  data,  as  seen  in  Figure  9. 

The  data  of  the  model  and  experiment  are  given: 

L  =  5/8"  D  =  1/8"  h  =  5 

D 


For  helium  gas 


J-‘  =  386 


ft 


V  =  5/3  . 


The  test  conditions  are  Moo  =11  T*  =  535°d  n  Ann 

At  R,  and  p  t  =  400  psia. 

Reynolds  number  per  inch  based  upon  its  conditions  at  the  edge  of  cone 
boundary  layer 

Rec /in  =  5  .  9  .  105 

Flow  parameters  at  the  edge  of  cone  boundary  layer  are  computed  by  us 
Reference  [31]  as  follows: 


Mc  =  6.46 


535 


1  +  'Ll.!. 
2 


M, 


u. 


=  /  7  8 

=  M 


=  35.  88°R 


=  144 


=  862  ft/sec 


c  •  ac 


5560  ft/sec 


T„  = 


Y 

(1  +•*" 1  mJ 

l 

-  Poo  Pc 

Ptc 

Ptc 

Pt® 

=  Tt 

1  +  v  ~  1 

C 

M2 

_  Px 

=  0. 

400.144 

11,000 


5.24 


535  =  12.93  R 

41. 33 


f  a  t0 


Pc  =  P®/  0.  35  =0.932.10"4 


Sc  = 


“c 

5.9.106  12 


7.85.  10'4  ft2 


=  fc-  Vc  =  7.31.10'®  1-^£ 

The  parameters  of  the  cavity  flow  are: 
Pr-av  =  Pe  =  O-85  •  Pc 


[Ref.  25] 


taking 


P«t  =  Pct 

-- --av  =  pcav  Pcavt  pc 

Pc  -  •  - L  •  — ct 

Pcavt  Pct  Pc 


-  (1  +  M  2  ) 


Me  =  6.7  Te  = 

L  nr  lrt-8  lb  sec 
Ve  =  6. 85.  10  — 

f  =  0.  826.  10"4  lb  sec' 


Te  =  33.8  R, 
•8  lb  sec 


}  e 


a  =  840  —  .  u.  =  5640  ft/sec 


sec  ) 


Assuming  that  the  reattachment  begins  at  x  =0.8  L,  [25],  and  at  this 


position  u  =  0.  587,  the  pressure  jratio  at  the  reattachment  point  is 


x  r 

Pcav 


1  +  1S-1  M2 

/  e 


1  +(i  -  u*2) 

z 


2.  65 


Pcav 


=  1.5 

0.85 


1.76 


From  Figure  4  of  Reference  [25],  [also  Figure  9  of  this  Report],  it  is 
estimated  that 


thus  ^  »  0.665 


This  result  shows  that  the  isentropic  compression  occurs  upstream  of 
reattachment. 

For  the  selected  model  configuration,  it  is  likely  that  8  <  ^  t*lus  the 
velocity  distribution  in  the  reattachment  zone  is  from  equation  (16). 

U(*  ^  =  \  ^  r o  +  fl  )  (L-x)  tan  8 

A  small  term  of  0  (L-x)  of  equations  (16)  is  neglected). 

Assuming  that  1  n  .n 

2  ( i  o  1  )  tan  8  =  C  i 

Then  ud  (x)  =Cl  (L  x) 

This  unknown  constant  Cj  is  determined  by  matching  with  the  experimental 
data  at  the  beginning  of  reattachment,  i.e.,  at  the  position  of  L-x  =0.2  L  =  | 

"95 


Cj  =3.17.10 


sec.  ft. 


Hence 


ud  =  3.  17.  105  (L  -  x) 


Since 


-i  =i  +i_i 

Te  2__ 

ad  =  144  .pTa 


m2  (  i  .  u  ) 


ad  =  144  \J  Ta  =  2750  ft/sec 
and  Md  =  U  d  /  ad 
Md  =  1.15.  102  <L  *  x) 


(21) 
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For  the  trajectory  of  F  -  ,  the  following  values  o'  F  and  <.  are  taken 

from  Reference  [5]  at  the  separation  and  reattachment  points,  namely, 

Fs  =2.85  k  s  =  0.63 

FR  =  .59  *R  m  0.693 

The  values  of  F  and  at  reattachment  point  are  those  of  Blasius  flow.  At 
the  point  of  beginning  of  reattachment. 

<o  -  1  - 

By  the  method  of  trial  and  error,  the  values  of  C  and  <  0  are  determined. 

Taking  C  =  15, 

0  1 . 25 

Z  *  “2  aiw*  x  =  0.8.L, 

it  becomes 


=  0.  83 


Figure  7  shows  the  completed  trajectory,  with  the  local  values  of  F  and 


along  the  reattachment  distance.  Since  local  values  of  F  and  K  are  known 

.  2(  u,.  1) 

as  shown  in  Table  1  the  local  values  of  A  (sJ J  p  (n  -  and  pressure  in  the 


reattachment  zone  are  computed. 


Table  1 

Reattachment  Pressure  Rise  (Computed) 


x 

L 

< 

F 

-  2  (uu  - 

f  (n 

(*) 

n 

1) 

-  u>) 

^Lm/ 

P 

p 

r  cav 

P 

Pc 

0.  2 

0.  83 

2.  85 

1 

0.  480 

1 

0.  85 

0.  16 

0.  802 

2.6 

1. 0255 

0. 3072 

1.0151 

0. 8628 

0.  12 

0.  775 

2.  35 

1. 0544 

0.  1728 

1. 0202 

0. 8671 

0.  1 

0.  765 

2.  23 

1. 0668 

0.  120 

1. 0180 

0. 8653 

-2. 


As  *een  in  Figure  9,  the  computed  values  of  pressure  are  in  excellent 
agreement  close  to  upstream  half  distance  in  the  reattachment  zone.  This  re¬ 
sult  confirms  that  C  =  15  is  an  acceptable  value  in  the  region  of  upstream 
reattachment  zone  of  a  laminar  cavity  flow  which  Glick  [5]  used  successfully 
for  the  region  of  reattachment  zone  downstream  of  shock  impingement. 

From  Figure  9  it  is  also  noticed  that  the  measured  value  of  pressure  reac 
close  to  its  cone  value  at  point  A  which  is  located  at  three -fourth  downstream 

distance  of  the  reattachment  zone.  This  value  of  pressure  corresponds  approx 
imately  to  that  computed  by  isentropic  compression,  with  the  Md  at  that  positic 

multiplied  by  the  average  value  of  ^  av  =  1  +  0.  665 _  =  0.8325 


between  separation  and  reattachment. 
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3.  Heat  Transfer  in  the  Reattachment  Zone 


At  present  no  analysis  for  the  heat  transfer  at  the  reattachment  point  is 
available.  However,  since  the  process  of  reattachment  is  similar  to  that  of 
stagnation,  for  the  evaluation  of  reattachment  heat  transfer,  the  known  stagna¬ 
tion  heat  transfer  formula  may  be  referred.  Therefore,  it  is  of  interest  to 
compare  the  empirical  data  of  reattachment  heat  transfer  with  those  predicted 
by  the  stagnation  heat  transfer  equation. 

At  hypersonic  speed,  Chung  and  Viegas  [15]  employed  the  following  semi 
empirical  formula  of  average  heat  transfer  in  a  reattachment  region  of 

0  <  —  <  1  based  upon  the  stagnation  heat  transfer  analysis  of  Lees 

L 

[33  ]  for  highly  cooled  wall,  namely 


Pf  ^  J  j  iT  ~  ^  (°»76  +  l.Ull  fe  ) 


(22) 


Pt 


This  equation  inapplicable  for  two-dimensional  flow  with  fj  =  45  and 


also  for  0.  1  — 


<  Pe  < 


Pt 


—  0.  5  where  H  is  total  enthalpy. 


If  this  equation  is  assumed  to  be  applicable  for  the  axial  symmetric  flow, 


—  =  -  =  0.57  and  for  small  area  of  reattachment  similar  to  Nicoll's 

pR  1.76 

experimental  model  [Figure  8b.  Ref.  [25  ]]  ,  since 

q  =  h (Tw  -  Tawj 


equation  (22)  may  be  written 


H. 


H 


*R  = 


w 


■  w 


-  T 


aw 


M  R 


0.  69 


assuming 


lUs-Is  fg  {£_  «  Pg  t,  TR 

**•.  *a  /» _ f.  Vc  „  ,,,Pa  ,,Tr  5 


Me  t  /  f  fle  Pe  Tr  u  P« 

and  at  the  reattachment  where  y  =  0.  665  taking  Tj^  s  Taw,  it  becomes 

“h  -  «•*  Jf.%  c„ 


"*  .  0.9US  y  ?, «.  J  ^  Cp 

For  cone  . 

he  *  8tc  *  /c  *  Uc  •  P 
and  at  the  reattachment  point, 

Ba0  -  5,9.1c)5  (1.25  +  |  )  -  1.108.106 

For  this  Reynolds  number,  from  Figure  2,  Reference  [25]  , 

St  s  5.10-1* 

an-j 

?c  uc  =  0.932  .  10"*  .  C 
Therefore  P 

hc  =  2.64.  10  .  Jcp. 

From  equation  (23)  and  (24).  For  Ujj  =  0.  587.  5640  *  330  ft/sec. 


0.  505 


Thus  hR/hc  is  inversely  proportional  to  frR. 

Chung  and  Viegas  formula  yields 

(a)  For  Lr  =  0.  008-  L  =  0.  008  |  =  0.00417  ft 


7.82 


(b)  ForLR  -0.3.  L  -  0. 2-  |  .  ^  ft. 


=  4.  95 


of  Case  (a)  corresponds  to  the  measured  value  of  Lr  (Ref,  [15  ]for 
~  45°  and  Lr  of  Case  (b)  is  that  of  the  measured  Lr^  =  0.  2»  L  by  Nicoll  [25], 


By  comparing  these  average  values  of  hR/hc  of  equation  (26),  with  Nicoll's 

(Figure  8b,  Ref.  [  25  ]  )  average  value  of  hR/hc  *  2.63. 

Case  (a)  (hR/hc)  Chung  / 

Viegas 

Case  (b)  (hR/hc)  Chung  /  /  hR  \  *  1.88  ^  2 

Viegas  —  ) 

\  hc  /  Nicoil 

Hence,  Chung-Viegas'  semi-empirical  formula  for  reattachment  heat  trans¬ 
fer  predicts  higher  heat  transfer  values  than  those  obtained  by  experiment. 

In  order  to  more  accurately  establish  the  prediction  of  heat  transfer  in  the 

reattachment  zone,  the  distribution  of  velocity  u__  along  the  reattachment  surface 
__  s 

s  and  the  velocity  gradient  at  the  reattachment  point  must  be  determined.  Because 

the  flow  process  at  the  reattachment  point  is  complex,  the  reattachment  velocity 

gradient  is  determined  semi-empirically,  but  u_  ,  the  potential  flow  velocity  along 
_  s 

s  downstream  of  reattachment  point,  is  computed,  using  Squire's  model  of  cavity 
flow  involving  a  core,  and  its  surrounding  boundary  layer  as  shown  in  Figures  1 
and  2. 

Although  the  viscosity  is  responsible  for  the  development  of  the  motion  in 

the  cavity,  within  the  core  the  viscosity  effect  is  small.  Thus,  the  flow  along  the 

outside  border  line  whirh  is  also  the  border  line  of  the  boundary  laryer  is  assumed 

to  be  inviscid.  Then,  u_  the  potential  .velocity  along  the  wall,  can  be  computed 

s 

along  this  border  line  of  the  boundary  layer. 

Hence,  by  simplicity  assuming  the  flow  motion  in  the  reattachment  zone  is  in- 

fc 

compressible,  the  following  Bernoulli's  equation  is  applicable  in  the  reattachment 


hR  \  *  2.975  3 

hc  )  Nicoil 


zone. 


Pr 

? 


T 


Then 


/ 


PR-P 

f  ud2, 


P8  -P 


Id 


!27) 


where 


"d  *  lv 

2 

The  heat  transfer  in  the  region  of  stagnation  point  of  cold  blunt 


nose  is  givf 


by  Lees  &2j.  Writing  this  equation  for  ud  along  the  dividing  stream  line  whic 
may  be  considered  the  upstream  velocity  of  reattachment,  similar  to  upstream 
free  stream  velocity  approaching  stagnation  point,  one  obtains 


u 


-  j_,(P/Pp)  ^  ?c 


c 

/ 


u 

P_  a  2  - 

%  ud  r°  ds. 


u. 


du_ 

s 

d  5 


'll  »  0 


(28) 


where  rQ  is  radius  of  cross  section  of  body  of  revolution.  Equation  (28)  is 


a  function  of  urf  because  the  heat  transfer 


process  in  the  reattachment  zone, 


similar  to  the  case  of  pressure 


rise  in  the  reattachment  zone  can  be  formulatec 


by  the  flow  properties  along  the  dividing  stream  line. 


Thus 


d  u 


d  .? 


constant 


the  constant  is  to  be  evaluated. 
Equation  (28)  reduces  to 

a  -Up/Q 


p0-p 


*  r. 


% 


[ 


B 


SLlf  r2 
qd  0  5 


1 

2 


(29) 


Hence  by  knowing  pressure  distribution  along  the  dividing  stream  line  with¬ 
in  the  reattachment  zone  for  a  given  geometry,  the  heat  transfer  in  the  reattach¬ 


ment  zone  can  be  computed  if  Cj  is  determined  empirically. 

Numerical  Evaluation  of  Heat  Transfer  in  the  Reattachment  Zone 


Chung  and  Viegas  Cl 5]  found  the  stream  function  by  solving  the  reattachment 
flow  as  incompressible  two-dimensional  inviscid  but  rotational  as 
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where  A  is  a  parameter  and  its  numerical  value  is  ^  =  1  Ra^ 

2.222  — 

u  L 

and  R.  =  ® 

T.  — - — 


Since  u 


H 


the  velocity  gradient  at  the^reattachment  point  is 


Now  with  sufficiently  large  n  value,  say  nj 


1*  (-2-)2  »2 


(-1— )*  n2 


hence  the  series  for  terms  larger  than  n^  become 


2  E 

“1 


2  \  L 


R  E 


_ s  - a? 

TT  w  -  o 

x  "fc  L 

This  singularity  of  velocity  gradient  at  the  reattachment  is  physically  imp' 

sible.  As  previously  mentioned  Burggraf's  [  13  ]  analysis  leads  to  the  reattach 

ment  heat  flux  to  square  root  singularity. 

At  present  no  other  analysis  of  reattachment  velocity  gradient  is  available 

therefore  its  value  is  evaluated  empirically  from  experimental  data  of  Nicoll 

[25  ]with  the  test  model  sketched  in  Figure  8. 

Since  with  this  model,  the  initial  point  of  pressure  and  heat  transfer  coinc 

approximately  at  1  -  —  =  0.  2,  it  is  simple  to  correlate  the  pressure  and  heat 

L 

transfer  in  the  reattachment  zone. 

Rewriting  equation  (29)  by  cone  values  of  pressure. 


P  P, 


f  1  !s  i/  If  f  -  L  \ 

J  pc  %  If  Id  v  Pc  pc  / 


The  heat  transfer  coefficient  ratio  is, 


i'o2  dT 


h  =  h 


St  (Tw  “  Taw) 


(31) 
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The  local  value  of  h^/hc  is  estimated  by  multiplying  a  factor  of  1.15  to 
the  average  value  of  h/hc  in  the  small  area  at  the  reattachment  point,  Figure  8b, 
Ref. [25 ]  . 


Thus,  hH  * 

T*  ~  1*5  (  J*  av  )  *  1.15  .2.63  3.0 


Since, 


T  ,  T 
d  e 


[  l  *  1?  V  (l-u **}]  -  36 


and  f.  «*  pd 
d 


U2  .  0.85 


g^Td  32.2  .  386  .  365 
ud  *  0.587.56U0  *  3300  ft/sec 


0.792.10“^  ib  88C2 


ft 


and 


<ld  =  Pd  “d  =  U3 


it  becomes  pc/qd  2b  1. 

X 

The  value  of  C_  is  evaluated  by  equation  (30)  at  a  position  of  1 - =0.  04, 

^  Li 


where  —  =  1.34  or  9W 


0.446. 


>wr 


The  values  of  denominator  and  nominator  of  equation  (30)  excluding  C2  are 
listed  in  Table  2. 


Table  2 

Computed  Values  of  Equation  30 


X 

L 

I  £.  h.  ,/ p=  -p  \ 

*  H  [pT  frj 

1 

f  ILfsi/—  /  Pr  - p  )  *„2<« 

-'o  p0  PRi<id  (_?;  id  0 

0.  2 

0.  3425.  10“2 

2.  7985. 10'6 

0.  16 

0. 3525. 10"2 

2.5335. 10"6 

0.  12 

0.  365.  10”2 

2.  326.  10’6 

0.  08 

0.  406.  10"2 

1.965.  lO-6 

0.  04 

0.  4515. 10“2 

1.54. 10'6 

0 


0 


0 


Since  1  -  =0*  04, 


it  follows  °2  *  (1.5U.10"6)*  = 


0.1*lj6 


=  l.QL.KT2 


or  ^  1 

Denoting  d  =  2D,  the  stagnation  velocity  gradient, 


i  -° 


|  i-  66.5  -  1.39 
8  12 


From  Figure  7-. 1 6,  Ref.  [33  ]  it  is  seen  that  this  value  of  stagnation  veloci 


gradient  corresponds  to  that  of  blunt-nosed  body  for  incompressible  flow  as  we 
as  to  Newtonian  flow  for  supersonic  flow. 

Finally,  the  local  values  of  h  is  determined  from  equation  (31).  Since  th 


ratio  of  local  value  of  (T»  ~  Taw)R 


is  not  known,  it  is  necessary  to  evali 


this  ratio  first. 


Tw  “  Taw 


At  reattachment,  1  -  —  =0,  (Ty  _  ^sw)r  =  1  and  at  the  initial  point  of 

L  Tw  -  Ta*  b 

reattachment,  i.  e. ,  at  1  -  ~  =0.2,  q  anc*  =  ^5 


reanacnment,  l.  e.  ,  at  l  -  —  -  a  or--.  v, 

L  =  o.251  nc 

(Figure  8b,  Reference  [  25 ]  ),  thus  from^quation  (31)  (Tw  -  Taw)^ 


=  0.4 


Assuming  a  linear  variation  of  ^w  ”  Taw)g 

Tw  —  Taw 


Tv  “  Taw 


in  the  reattachment 


distance,  thus  drawing  a  straight  line  which  has  its  initial  and  end  values  of 


Tw  -  Taw)R 

Tw  "  Taw 
( 1  w  ”  Taw)R 
Tw  *  Taw 


equal  0.465  and  1,  as  seen  in  Figure  10,  the  local  values 


are  determined  as  shown  in  Table  3  or  Figure 


Then  from  equation  (31)  local  values  of 


are  computed  as  indicated  in 


Table  3  and  Figure  10. 


-32- 


Tablt  3 

Computed  Values  of  Equation  (31) 


1  -  - 

.  L . 

<T“  '  I")*/  (Tw  -  *„) 

hAc 

0.  2 

0.251 

0.465 

0.35 

0.  16 

0.271 

0.  57 
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0.  12 

0.292 

0.  676 

0.59 

0.08 

0.354 

0.  782 

0.828 

0.04 

0.446 

0.89 

1.  19 

0 

1 

1 

3 

The  comparison  of  computed  values  of  h/hc  with  those  measured  shows  good 
agreement  among  them,  as  seen  in  Figure  10.  Therefore,  it  appears  that  the 
assumed  condition  of  incompressible  flow  prevails  in  the  reattachment  zone. 


9? 

In  order  to  confirm  this,  M  and  density  ratio  of  — —  are  computed  with 

?  f’d 

the  temperature  distribution  shown  in  Table  3. 

Since 


u  u 

J  =  JS  H 


ud 


and 


W 


Md  *  1.212,  with  teat  condition  M_  =  0.812  -  0.98  <  1 

a 

in  the  reattachment  zone. 


Then 


{ \ 


*JL  3s  = 

pd  * 


1.11 


in  the  range  near  80%  of  the  reattachment  zone  extending  from  the  initial  point 
of  reattachment  to  downstream,  distance. 

These  figures  of  M_  and  p~  may  justify  the  assumed  condition  of  in- 

s 

compressible  flow  in  the  reattachment  zone. 
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Assuming  that  the  known  stagnation  heat  transfer  equation  for  blunt  body 
is  applicable  for  the  reattachment,  the  following  formula  [34]  is  used. 

\l  =  f2|/|"  jyJ  /d  J*d  CP  (Tw  “  Taw^ 

8=0 


where  =  0.  57.  P f  for  two-dimensional  flow 

-0.6 

f 2  =  0.  763.  Py  for  axially  symmetric  flow 
The  equation  (33)  is  written  by  flow  property  along  the  dividing  stream  line 
at  the  beginning  of  reattachment,  i.  e. ,  £d  ^  (For  stagnation  heat  transfer 
for  blunt  nose  exposed  free  stream,  this  term  is  f  ^  ). 

Since  f>  ^  ^  B^e  and  if  one  assumes  a  linear  relation  of  -  -p  equation 
(33)  becomes 


/[5j,=0 


This  equation  expresses  hj^  based  upon  the  flow  condition  at  the  outer  edge 
of  viscous  layer. 

Because  hc  =  Stc.pc.uc.c  as  before, 


JL  V^Me 

St«  fc  » 


For  the  test  condition  of  the  model  sketched  in  Figure  8, 

l  ,  1 

(  >  =  (0.826.10"^  .  6.85.10’8)  1  =  2.377.10 


377.10"6  ClAee^)  2 


fc  uc  =  0. 932. 10'*  .  5560  = 
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Stc  *  5.  10  referring  to  the  cone  value  at  the  reattachment  point,  and 


sine* 


77  (— lL  *  66.5 

“ <*  v  «»  A --o  »  Ai  /  i 


-  22. 104 


from  equation  (35) 

A  *  f2  ^»377«10‘’6  ^  22.10^  *  f2  .  4. 3 

hc  5.10^  .  0.52 

In  helium  Pr”^*  ^  =  (0.  68)  ^  =  1.26 


f2  ***  0.72 

for 

two-dimensional  flow 

f2  0.96 

for  axially  symmetric  flow 

Finally 

hR 

hc 

3.  12 

for  two-dimensional  flow 

hR 

^c 

4.  15 

for  axially  symmetric  flow 

but  hj^ 

hc 

3.0 

from  experiment. 

From  these  results,  it  appears  that  reattachment  heat  transfer  may  be 
predicted  in  a  good  agreement  with  experimental  data  by  using  stagnation  heat 
transfer  equation  for  blunt  ncsed  body  assuming  reattachment  heat  transfer  is 
essentially  a  two-dimensional  process,  because  the  assumption  of  axially  sym 
metric  flow  could  lead  to  higher  value  of  reattachment  heat  transfer. 

Actual  reattachment  process  is  very  complicated  involving  geometry  of 
reattachment  surface,  angle  0  ,  reverse  flow,  etc.  Hence  further  study  is 

needed  to  clarify  the  difference  between  the  reattaching  and  stagnating  flow 

t 

processes  involving  the  velocity  gradient  parameter  at  the  reattachment  point. 
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4.  Discussions 

The  distributions  of  pressure  and  heat  transfer  of  cavity  flow  are  affected 
by  the  geometry  of  cavity;  i.  e. ,  by  L/D,  the  ratio  of  length  to  depth.  As  shown 
in  Figure  11  and  12,  the  pressure  rises  immediately  after  separation  and  contin¬ 
ues  to  rise  in  the  whole  region  of  the  cavity  if  L/D  is  large,  say  L/D  =  10,  but  the 
heat  transfer  rates  remain  constant  except  in  the  narrow  region  of  reattachment 
for  L/D  =  40/3.  On  the  other  hand,  as  shown  also  in  Figure  11  and  13,  if  L/D  is 
smaller,  say  L/D  =  5,  the  pressure  remains  constant  upstream  of  reattachment, 
only  rising  sharply  in  the  reattachment  zone,  but  at  the  same  L/D  the  heat  transfer 
rate  reaches  minimum  at  the  center  of  cavity  and  rises  rapidly  approaching  the 
reattachment  point.  These  behaviors  of  pressure  and  heat  transfer  appear  strongly 
influenced  by  the  recirculating  flow  characteristics  depending  upon  the  depth  of  its 
cavity.  Because  the  initial  points  of  rise  of  pressure  nnd  heat  transfer  must  be 
known  for  the  prediction  of  reattachment  process,  the  location  of  the  initial  point 
of  reattachment  must  be  determined  for  a  given  L/D. 

Burggraf's  [  13  ]  paper  which  considers  the  various  depths  of  the  cavity  was  ob 
tainedat  a  late  stage  in  the  preparation  of  this  report.  Thus,  no  investigation  of 
these  reattachment  behaviors  has  been  made  using  Burggraf's  analysis. 

It  is  found  that  Glick's  method  [5]  for  computation  of  reattachment  flow  in- 

* 

volving  shock  interaction  is  applicable  for  the  cavity  flow.  Also  the  semi-empirical 
average  mixing  rate  correlation  parameter  C  is  the  same  for  these  two  flows.  But 
it  is  still  unknown  whether  Glick's  analysis  is  also  applicable  to  other  types  of 
separated  and  reattached  flow 3  such  as  caused  by  forward  as  well  as  reward  facing 


steps,  corners  and  ramps,  etc. 
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An  analytical  study  is  desirable  in  order  to  evaluate  the  reattachment  velocity 
radient  parameter  for  arbitrary  geometry  of  reattachment  surface,  angle  3 
etween  dividing  stream  line  and  reattachment  wall  and  reverse  flow  etc. 


5.  Conclusions 


An  investigation  of  the  cavity  flow  based  on  Chapman-Squire  type  physical 
model  referring  to  the  experimental  data  of  a  particular  limited  model  at  free 
stream  Mach  number  11,  leads  to  the  following  conclusions. 

Predictions  of  pressure  rise  and  heat  transfer  in  the  reattachment  zone  can 
be  made  by  flow  conditions  on  the  dividing  stream  line. 

Pressure  rise  can  be  computed  by  Glick's  procedure  of  analysis  with  an 
average  mixing  rate  correlation  function  C  =  15  in  the  mixing  zone  upstream  of  re¬ 
attachment  employing  Crocco-Lees'  mixing  theory.  Up  to  upstream  half  zone  of 
reattachment,  the  pressure  rise  takes  place  by  isentropic  compression,  but  ap¬ 
proaching  reattachment  point,  due  to  compressibility,  viscosity,  geometry  of 
reattachment  surface,  reverse  flow,  etc.  the  pressure  rise  is  less  than  that  of 
isentropic  compression. 

An  estimation  of  reattachment  heat  transfer  may  be  made  by  using  the  known 
stagnation  heat  transfer  equation  for  blunt  nosed  body. 

However,  additional  experimental  data  of  pressure  and  heat  transfer  distri¬ 
bution  in  the  reattachment  zone  of  various  geometrical  configurations,  as  well  as 
the  analytical  prediction  of  reattachment  heat  transfer  are  needed  to  clarify 
completely  the  reattachment  heat  transfer  in  the  cavity  zone. 
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Figure  7  F(K)  Trajectory  for  Laminar  Cavity  Flow 
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